To APPEAR IN THE ApJ 

Preprint typeset using IATgX style emulateapj v. 08/22/09 



(N 

o 

(N 
>Y 



in 

(N 

< 
O 

9 L,: 
6 

CO 

o3 



> 
O 

in 
i> 
in 

in 
o 

(N 



13 



FROM DUSTY FILAMENTS TO MASSIVE STARS: THE CASE OF NGC 7538 S 

Raul Naranjo-Romero 1 , Luis A. Zapata', Enrique Vazquez-Semadeni 1 , Satoko Takahashi 2 , Aina Palau 3 , and Peter Schilke 5 

To appear in the ApJ 

ABSTRACT 

We report on high-sensitivity and high-angular resolution archival Submillimeter Array (SMA) observations 
of the large ( ~15000 AU) putative circumstellar disk associated with the O-type protostar NGC 7538 S. Obser- 
vations of the continuum resolve this putative circumstellar disk into five compact sources, with sizes ~ 3000 
AU and masses ~ 1OM . This confirm the results of recent millimeter observations made with CARMA/BIMA 
towards this object. However, we find that from most of these compact sources eject collimated bipolar 
outflows, revealed by our silicon monoxide (SiO 7=5-4) observations and confirm that these sources have a 
(proto)stellar nature. All outflows are perpendicular to the large and rotating dusty structure. We propose 
therefore that, rather than being a single massive circumstellar disk, NGC 7538 S could be instead a large and 
massive contracting or rotating filament that is fragmenting at scales of 0.1 to 0.01 pc to form several B-type 
stars, via the standard process involving outflows and disks. As in recent high spatial resolution studies of 
dusty filaments, our observations also suggest that thermal pressure does not seem to be sufficient to support 
the filament, so that either additional support needs to be invoked, or else the filament must be in the process of 
collapsing. An SPH numerical simulation of the formation of a molecular cloud by converging warm neutral 
medium flows produces contracting filaments whose dimensions and spacings between the stars forming within 
them, as well as their column densities, strongly resemble those observed in the filament reported here. 
Subject headings: stars: formation — ISM: jets and outflows — ISM: individual objects (NGC 7538 S) 



1. INTRODUCTION 

The observational investigation of the formation of young 
and very massive stars (with L* > 10 s M , i.e. O-type) has 
been inherently a difficult task. These young stars are rare 
and found far away from us, at typical distances of more than 
2 kpc. Their evolutionary time scales are also much shorter 
than those found in low and intermediate mass stars (McKee 
& Tan 2002). At present, there is only a handful of very mas- 
sive star forming regions that have been extensively studied 
at infrared, (sub)millimeter, and/or radio wavelengths using 
sensitive, high angular resolution observations. Some of them 
include G240.3 1+0.07 (Qiu et al. 2009), IRAS 16547-4247 
(Garay et al. 2007; Franco-Hernandez et al. 2009), W5 1 North 
(Zapata et al. 2009a, 2010), NGC 7538 IRS1 (Qiu et al. 201 1), 
and G24 Al (Beltran et al. 2011). These young O-type stars 
are found in association with large (~ 5000 AU) accreting ro- 
tating structures that seem to energize powerful bipolar out- 
flows suggesting that these stars may form in a similar way to 
the low mass stars, i.e. via outflows and disks. However, it 
is very likely that those large accreting structures are forming 
more than one single high-mass star at once. These structures 
then are not disks in the sense of accretion disks that funnel 
gas to a single star. 

Here, we present high angular resolution archival obser- 
vations using the Submillimeter Array at millimeter wave- 
lengths towards the extremely large accreting structure sur- 
rounding the O-type protostar NGC 7538 S. This is a cold core 
with a size of 0.5 pc, and a dust mass of ~ 2000 M associated 
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with H2O, OH and CH3OH (but not SiO) maser emission, and 
weak free-free emission (Sandell & Sievers 2004; Pestalozzi 
et al. 2006; Zapata et al. 2009b). This region is at a distance 
of 2.65 + 0.12 kpc (Moscadelli et al. 2009). In the middle 
of this cold core resides a large (~ 15 000 AU) and massive 
(~ 400 M Q ) structure interpreted as a nearly edge-on rotating 
disk that surrounds a protostar with a mass between 20 - 40 
M Q (Sandell et al. 2003; Sandell & Wright 2010). However, 
the values reported by Sandell & Wright (2010), which should 
be the most recent estimates, are significantly lower, 85-115 
M for the putative disk, and 14-24 M for the enclosed mass 
of the central protostar and the inner, rotationally supported 
part of the disk. 

A single, young, and powerful outflow perpendicular to the 
rotating disk has also been mapped in molecular lines that 
includes SiO (7=2-1) and HCO + (7=1-0) (Sandell & Wright 
2010). The bipolar outflow is compact (0.2 pc), has a low 
velocity (Av ~ 40 km s~'), and has its blueshifted side towards 
the northwest, and the redshifted one to the southeast. 

Very recent BIMA and CARMA observations have re- 
solved the large and massive disk into four compact millime- 
ter sources (S aia2 ,b,c) (Wright et al. 2012; Corder 2008). 
Whether the massive and unstable disk fragmented in four 
pieces or whether a dusty core fragmented into four protostars 
is unclear. 

2. OBSERVATIONS 

The observations were obtained from the SMA 5 archive, 
and were made on October 2005 and October 2006, when the 
array was respectively in its compact and extended configu- 
rations. The independent baselines in the compact and ex- 
tended configurations ranged in projected length from 5 to 

5 The Submillimeter Array (SMA) is a joint project between the Smithso- 
nian Astrophysical Observatory and the Academia Sinica Institute of Astron- 
omy and Astrophysics, and is funded by the Smithsonian Institution and the 
Academia Sinica. 
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TABLE 1 

Observational parameters of the SMA compact sources 





Position" 


Flux Density 


Flux Peak 


Deconvolved Angular Size'' 


dust-gas mass c 


Name rf 


ff(J20()0) 


<5(J2O00) 


(mJy) 


(mJy) 




(Mo) 


SMA1 (CARMA S c ) 


23 13 44.483 


+61 26 47.70 


130±30 


41±4 


i:'5±0:'6x i:'2±0"8; 78° ±40° 


10.0 


SMA2 (CARMA S„) 


23 13 44.770 


+61 26 48.85 


110±26 


44±4 


i:'3 ± o:'4 x o:'4 ± o:'5; 19° ± 53° 


8.0 


SMA3 (CARMA Sai) 


23 13 44.979 


+61 26 49.41 


143±34 


56±3 


i:'5±o:'5 x i:'3±0"6; 45° ±35° 


13.0 


SMA4 (CARMA Sm) 


23 13 45.099 


+61 26 50.04 










SMA5- 


23 13 45.267 


+61 26 50.47 


80±30 


12±3 


i:'6±o:'9x i:'o±o:'8; 133° ±44° 


7.0 



"Units of right ascension are hours, minutes, and seconds and units of dec- 
lination are degrees, arcminutes, and arcseconds. 

b Major axis X minor axis; position angle of major axis. The values were 
obtained using the task JMFIT of AIPS. 

c The values of the masses obtained here have uncertainties of order a fac- 
tor of 2 due to the error in the estimation of temperatures of the millimeter 
sources and the error in the dust mass opacity coefficient at this wavelength. 
See main text for dust temperature and dust opacity assumptions. 

d The CARMA names were obtained from Wright et al. (2012); Corder 
(2008). 
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Fig. 1 . — Color and grey contour map of the 1 .3 mm continuum emission from NGC 
7538 S region obtained with the extended configuration and overlaid in black contours 
with the millimeter continuum image obtained with the compact configuration. The grey 
contours range from 20% to 90% with steps of 10% of the peak of the emission. The 
black contours range from 55% to 90% with steps of 7% of the peak of the emission. 
The synthesized beams of the 1.3 mm continuum observations are described in the text 
and are shown in the bottom left corner of the image. 



50 kA and 15 to 120 kA, respectively. The phase reference 
center for the observations was at R.A. = 23hl3m44.46s, 
decl.= +61°26'49.8" (J2000.0). The frequency was centered 
at 217.238539 GHz in the Lower Sideband (LSB), while the 
Upper Sideband (USB) was centered at 227.238539 GHz. 
The primary beam of the SMA in the neighborhood of 220 
GHz has a FWHM ~ 50". The emission from the whole NGC 
7538 S region falls well inside of the FWHM. The largest 
brightness distribution to which these observations are sensi- 
tive enough (more than the 50% level) is ~ 27" in FWHM 
when the distribution is Gaussian (Wilner & Welch 1994). 

The SMA digital correlator was configured in 24 spectral 
windows ("chunks") of 104 MHz each, with 128 channels dis- 
tributed over each spectral window, providing a resolution of 
0.812 MHz (~ 1.106 km s 1 ) per channel. 



The zenith opacity (t23qg#z)> measured with the NRAO 
tipping radiometer located at the Caltech Submillimeter Ob- 
servatory, was between 0.1 and 0.2, indicating reasonable 
weather conditions during the observations. Observations of 
Uranus provided the absolute scale for the flux density cali- 
bration. The gain calibrators were the quasars BL LAC and 
0102+584. The uncertainty in the flux scale is estimated to 
be between 15 and 20%, based on the SMA monitoring of 
quasars. Further technical descriptions of the SMA and its 
calibration schemes can be found in Ho et al. (2004). 

The data were calibrated using the IDL superset MIR, orig- 
inally developed for the Owens Valley Radio Observatory 
(OVRO, Scoville et al. 1993) and adapted for the SMA. 6 
The calibrated data were imaged and analyzed in the standard 
manner using the MIRIAD and KARMA (Gooch 1996) soft- 
wares. We set the ROBUST parameter of the task INVERT 
to to obtain an optimal compromise between resolution and 
sensitivity. The resulting rms noise for the continuum high 
resolution image was around 10 mJy beam -1 at an angular 
resolution of L'33 x L'04 with a PA. = -63.8°, while for 
the low resolution image it was around 5 mJy beam -1 at an 
angular resolution of 3"23 x 3" 15 with a P.A. = 6.0°. 

For the line emission, the resulting synthesized beam was 
L'65 x L'37 with a P.A. = -69.8° (SiO). The rms noise for the 
line image was about 150 mJy. We concatenated the uv-data 
from both configurations to make the line emission map. It 
is important to mention that in these observations many more 
lines were detected. However, in this work we concentrate 
only on the molecular distribution of the SiO. This molecule 
is well known to be good tracer of high density gas in outflows 
with critical densities of more than 10 3 cm -3 . 

3. RESULTS 

3.1. Continuum emission 

In Figure 1, we show the resulting millimeter (A = 1300 
/mi) continuum image obtained with the SMA in the compact 
and extended configurations from the high-mass star form- 
ing region NGC 7538 S. From the compact configuration, we 
detect a single northeast-southwest extended dusty source re- 
ported for the first time by Sandell et al. (2003); Sandell & 
Wright (2010) and interpreted as an extremely large edge- 
on circumstellar disk. However, Wright et al. (2012) made 

6 The MIR-IDL cookbook by C. Qi can be found at http://cfa- 
www.harvard.edu/~cqi/mircook.html 
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Fig. 2. — Integrated intensity color and contour map of the SiO (7=5-4) line emission observed toward NGC 7538 S overlaid in contours and green color with the SMA continuum 
emission from Figure 1. The blue and red contours are from 30% to 90% with steps of 7% of the peak of the line emission; the peak the SiO(/=5-4) emission is 100 Jy Beam -1 km s _1 . 
The white contours are from 30% to 90% with steps of 7% of the peak of the emission; the peak of the 1 .3 mm map is 730 mJy BeanT ] . The green contours are from 20% to 10% with 
steps of 7% of the peak of the emission. The synthesized beams of the SiO and 1 .3 mm continuum images are described in the text and are shown in the left lower corner of the figure. 
The white dot marks the position of the infrared source found in the archive of the Spitzer Space Telescope. The white arrows mark the orientation of the collimated outflows revealed 
in the SiO (J=5-4) observations. 



BIMA and CARMA observations and resolved the large and 
massive disk into four compact millimeter sources (Sai,a2,b,c) 
(Wright et al. 2012). The nature of these source is still in 
debate. 

This source has a deconvolved size of 9"4+0"6x5"l +0"3 
with a position angle 58° + 8°. Its total flux density at this 
wavelength is 4. 1 + 0.5 Jy. On the other hand, in the extended 
configuration, we resolved the elongated dusty source into 
five compact sources (SMA1-5), with the compact objects be- 
ing located along the major axis of the source detected with 
the compact configuration. These compact sources have sizes 
of about 1.5" (or 3000 AU at the distance of NGC 3875 S, see 
Table 1). The deconvolved parameters of the sources SMA3/4 
were determined by fitting a two-dimensional Gaussian func- 
tion to the observed flux distribution using the task JMFIT 
of ALPS 7 . SMA3 was well resolved using this method while 
for the fainter source SMA4 failed. Thus, we only give the 
position on the sky for SMA4. A better estimation of the de- 
convolved sizes of SMA3/4 (CARMA S^i and S^2) are given 
in Wright etal. (2012). 

7 Astronomical Image Processing System (AIPS) of the NRAO: 
http://www.aips.nrao.edu/index.shtml 



Four of these sources (SMA 1,2,3 and 4) coincide well with 
the positions of the four compact sources found by Corder 
(2008); Sandell & Wright (2010); Wright et al. (2012) using 
observations of CARMA/BIMA, which they propose to cer- 
tainly have a protostellar nature. The source SMA 5 does not 
have any counterpart at mm wavelengths as revealed by the 
Figure 8 of Wright et al. (2012). We are thus caution about 
if this source is real or is an artifact by poor the uv-sampling. 
However, as our synthesized beam is larger to the ones ob- 
tained by the CARMA and BIMA observations in Wright et 
al. (2012), perhaps, we are only recovering some more ex- 
tended emission and revealing the non-compact sources. The 
density fluxes for the sources SMA1, 2, 3 and 4 that we ob- 
tained here were larger (by a factor of more than 3) than those 
obtained by Wright et al. (2012), again probably because we 
are recovering more extended emission. The deconvolved 
sizes are in good agreement with those presented in Wright 
etal. (2012). 

Following Zapata et al. (2012), we adopted a value of /(1.3mm 
= 0.1 cm 2 g -1 (Ossenkopf & Henning 1994) and assumed op- 
tically thin, isothermal dust emission and a gas-to-dust ratio of 
100, with a dust temperature of 30 K (Sandell & Wright 2010) 
for these millimeter compact objects, deriving dust masses 
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of about 10 solar masses (see Table 1). From our compact 
configuration observations, we estimate that the elongated 
northeast-southwest dusty structure reported by Sandell et al. 
(2003) and Sandell & Wright (2010) has a mass of 350 M Q , 
which is in very good agreement with the mass obtained by 
those authors. The assumption of a dust temperature of 30 K 
appears justified, see, e.g., Sandell & Sievers (2004); Sandell 
& Wright (2010). For example, Zheng et al. (2001) find am- 
monia to be optically thick in the direction of NGC 7538 S 
with a temperature of 25 K, however, the low resolution of 
the observations did not allow to measure the temperatures of 
every single mm component of NGC 7538 S. 

We noted that Wright et al. (2012) obtained a mass of 60 
M for the source CARMA S A i or SMA3. This value is very 
different from that we obtained (about 10 M ). This discrep- 
ancy could arise from the low value of /c 3mm = 0.002 cm 2 g" 1 
they used. There is not a general consensus of using a singu- 
lar value, however, this dust mass opacity at this wavelength 
seems to be too low.. 

3.2. Line emission 

In our 2 GHz lower side band of the observations, we de- 
tected the line SiO (.7=5-4) at a rest frequency of 217.104 
GHz. In Figure 2, we present a map of the intensity, inte- 
grated over velocity (moment 0), of the SiO line emission ob- 
served towards NGC7538S. The map is additionally overlaid 
with the two maps of the continuum emission obtained from 
the compact and extended configurations (Figure 1). The ve- 
locity integration is over the velocity ranges: -70 to -57 km 
sT l (blue) and -53 to -40 km sT l (red). This map reveals 
at least two collimated and bipolar outflows emanating from 
NGC 7538 S, all with similar orientations. The bipolar out- 
flows emanate with their blueshifted sides toward the north- 
west while the redshifted sides to the southeast. These out- 
flows have position angles of approximately -45° and with 
collimation angles of less than 10°. The extension of the out- 
flows is about 0. 1 pc. 

There is one more bipolar outflow, revealed here for the 
first time, which is located about 7" north of NGC 7538 S 
(Figure 2). However, we do not find any millimeter contin- 
uum source in the outflow center. The outflow is also very 
collimated and with a similar orientation to those found in 
NGC 7538 S. An infrared study with the archive images of 
the Spitzer Infrared Telescope revealed an infrared (3.6 ^m) 
faint source in the middle of the outflow, whose position is 
RA. = 23hl3m45.121s,decl.=+61°26'57.01" (J2000.0). We 
suggest that this infrared source could be its exciting source. 
For the rest of the outflows in the field, the compact millimeter 
sources reported here (SMA1, 2 and 3) appear to be exciting 
them. We also find an IRAC (4.5 and 5.8 /im) source asso- 
ciated with the millimeter compact source SMA3 & 4. The 
infrared emission detected towards this object might suggest 
a more evolved evolutionary state compared with the rest of 
the compact millimeter sources. However, the nature of this 
infrared source is uncertain and an in-deep infrared study of 
the SMA3 & 4 is beyond of the scope of the present Letter. 

The collection of collimated outflows reported here appear 
to be part of the massive, bipolar, unresolved outflow reported 
by Sandell et al. (2003) and Sandell & Wright (2010). This 
has similar position angle (+148°), range of radial velocities 
(-80, -40 km s _1 ), and projected size (25" or 0.2 pc). 

The compact source SMA3/4 is additionally associated 
with a radio thermal jet with an orientation similar to the 
molecular outflows (Sandell et al. 2003; Sandell & Wright 



2010; Wright et al. 2012) also suggesting that this might be 
the exciting source of another outflow. 

Following Zapata et al. (2009a), and assuming local ther- 
modynamic equilibrium (LTE), that the molecular emission 
is optically thin, an excitation temperature of 50 K, and an 
abundance ratio of SiO/H2 equal to 1 x 10~ 7 , we can estimate 
the mass of the outflows for the SiO molecule in the transition 
7 = 5-4. We obtain a mass of about 10 M for each outflow. It 
is important to mention that the fractional abundance between 
the silicon monoxide and the molecular hydrogen varies from 
one star-forming region (SFR) to another, and thus the mass 
of the outflow can be over/underestimated. The value for the 
fractional abundance assumed here seems to be consistent for 
a few SFRs (see Ziurys & Friberg 1987; Zhang et al. 1995). 

For a mechanical force of F M = 10 M e 10 km s _1 /10,000 yr 
= 0.01 M Q km s _1 yr -1 , and from the correlation presented in 
Wu et al. (2004) for the outflow mechanical force versus the 
bolometric luminosity of the exciting source, we very roughly 
estimate a luminosity for the central powering sources of the 
order of 10 3 L , which corresponds to a massive B-type pro- 
tostar. This spectral type for the central star is in good agree- 
ment with that obtained from the dust emission, as we esti- 
mated on the last section. 

4. NUMERICAL SIMULATION 

Figure 4 shows a filament forming in a new numerical sim- 
ulation of molecular cloud formation by converging warm 
neutral medium (WNM) flows, similar to the one presented 
by Vazquez-Semadeni et al. (2007), but with 8 times higher 
mass resolution. The simulation will be described in detail in 
a separate contribution, but here we give the basic details. It 
was performed with the SPH code Gadget 2 (Springel 2005), 
with a prescription for "sink" particle creation first applied by 
Jappsen et al. (2005). According to this prescription, when the 
local density n at any location within the simulation exceeds 
a certain threshold density n tri r, the SPH gas particles within 
that region are replaced by a single sink particle, whose mass 
and momentum equal the sums of the masses and momenta 
of the removed gas particles. After being created, the sink 
particles behave as point masses, and can continue to accrete 
gas from a region of radius r acc around them, if such gas is 
gravitationally bound to the sinks. 

The present simulation contains 296 3 ~ 2.6xl0 7 SPH parti- 
cles, and starts at a volume density n — 1 cm -3 and T = 5000 
K, in a numerical box of 256 pc per side. The total mass 
in the box is thus 5.23 x 1O 5 M . In this box, we set up as 
initial velocity conditions two cylindrical gas streams of ra- 
dius 32 pc and length 1 12 pc, about to collide head-on at the 
plane x = 128 pc (i.e., at the middle of the box), with veloc- 
ity 9.2 km s _1 (for more details on this setup, see Vazquez- 
Semadeni et al. 2007). For the sink particle parameters, we 
choose nthr = 3.2 x 10 7 crrT 3 and r acc = 0.04 pc. 

The left panel of Fig. 4 shows the column density at time 
27.3 Myr in a 9 X 9 pc region on the (y,z) plane, near the 
center of the simulation, with a density integration over the 
range 1 18 < x < 130 pc. In this region, a long filament has 
formed by gravitational contraction, with several star-forming 
sites along it. 

The square in the left panel of Fig. 4 delimits the region 
shown in the right panel. This panel shows the smaller-scale 
filament arising within the larger one, immediately after form- 
ing three sink (stellar) particles, with masses ~ 3M each. 

The similarity between the filament+sinks system in the 
simulation and NGC7538 S is striking. The length spanned by 
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the three sinks, as seen in the right panel of Fig. 3, is ~ 0.1 pc, 
while the corresponding length in NGC 7538 S is ~ 6", corre- 
sponding to 0.08 pc. The sinks in the simulation, as well as the 
individual sources in NGC 7538 S are almost evenly spaced. 
The width of the filament containing the sinks in the simula- 
tion is also of approximately the same size, as is the width of 
the dust-emitting region in the source (Fig. 1). We can also 
estimate the column density of NGC 7538 S, and compare it 
to that of the filament in the simulation. The estimated mass 
of the region is 35OM . Assuming that it has the same exten- 
sion along the line of sight as it is wide on the plane of the sky 
(~ 0.08 pc), and a cylindrical shape, with length ~ 9" ~ 0.12 
pc, the resulting column density is N ~ 10 24 cm 2 , which 
compares very well with the range of column densities for the 
filament in the simulation (right panel of Fig. 3). 

Finally, it is worth remarking that NGC 7538 S is located 
within a larger filament spanning at least 3' (corresponding to 
~ 2.4 pc; Sandell & Wright 2010), while the filament from the 
simulation shown in the right panel of Fig. 3 can be seen in 
the left panel of the same figure to be part of a larger filament, 
spanning also ~ 2 pc (the filament containing the box in this 
figure). 

Note, however, that, in the simulation, this filament is in 
a state of gravitational contraction, both accreting mass in 
the direction perpendicular to its axis as well as contracting 
along this axis. Nevertheless, the simulated filament is seen 
to strongly resemble the observed one in both size scale and 
spacing of the newly-formed stars. 

5. DISCUSSION AND CONCLUSIONS 

The dimensions (~ 3000 AU), masses (of about 10 M Q ), and 
the fact that from most of the sources molecular collimated 
outflows or thermal jets emerge, suggest that those compact 
sources are probably large envelopes or perhaps dusty large 
disks associated with massive stars in the process of forma- 
tion. We thus conclude that NGC 7538 S is not likely to 
be a single, young O-type star surrounded by an extremely 
large disk, but rather maybe a large contracting or rotating 
filament that is fragmenting at scales of 0.1 pc to 0.01 pc 



to form multiple high-mass stars that involve outflows, en- 
velopes and disks. Similar molecular and dusty fragment- 
ing filaments have already been reported in the literature in 
OMC3 and G28. 34+0.06 (Takahashi et al. 2006; Wang et al. 
2011). For the case of G28.34+0.06, observations also from 
the SMA revealed a string of five dust cores of along the 1 
pc IR-dark filament called G28.34+0.06. The cores are well 
aligned at a position angle of +48°, and are regularly spaced, 
with an average projected separation of 0.16 pc. For OMC3, 
there is a large filament that also fragments at scales of 0.01 
pc to form the the protostellar object MM7 (Takahashi et al. 
2006). This large filament also shows large velocity gradients 
along it. This case is very similar to the case of NGC 7538 S 
reported here. 

Following Pillai et al. (201 1) and using the density derived 
from the emission detected in the compact configuration (1.3 
X 10 7 cm" 3 ), a velocity dispersion with a thermal line width 
(assuming 30 K) of 0.77 km s _1 (dispersion of 0.33 s _1 ), we 
estimate a Jeans mass of Mj = 0.76 M Q and a Jeans length of 
As = 1930 AU. Thus, for our case, the Jeans mass is smaller 
than the measured masses (by one order of magnitude, as in 
G28), and the Jeans length is slightly smaller than the mea- 
sured separation between the fragments. So we also conclude, 
in agreement with Wang et al. (2008, 201 1), that thermal pres- 
sure does not seem to be enough to support the source, which 
appears to be a filament, so that either additional support 
needs to be invoked, or else that the filament is in the pro- 
cess of gravitational contraction, as is the case of the filament 
in the numerical simulation presented here. 

Finally, we note that our results are consistent with re- 
cent BIMA, CARMA and SPITZER observations presented 
in Wright et al. (2012) and confirm the presence of multiple 
compact source located along of NGC7538 S. 
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